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Abstract 

The low-energy scattering of two charged particles is analyzed using a renormalization group 
approach based on dimensional regularization with power-divergence subtraction. A nontrivial 
solution with a marginally unstable direction is found, corresponding to a system with a bound 
state at zero energy. For purely energy-dependent perturbations around this solution, the power 
counting agrees with that from Wilsonian methods. These terms in the effective potential are in 
direct correspondence with the the terms in the Coulomb-distorted effective-range expansion. We 
also study perturbations that depend on off-shell momenta as well as energy, and we show that 
these affect only the off-shell form of the scattering matrix. These terms are of higher order that 
the corresponding energy-dependent ones and so terms in the potential that depend only on the 
off-shell momenta do not have definite orders in power counting. 
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I. INTRODUCTION 



The effective field theories (EFTs) that have been developed in recent years are now 
important tools for analyzing the scattering of low-energyparticles, particularly in the con- 
text of nuclear forces. (For reviews, see, e.g., Refs. Q, [s].) They can provide systematic 
expansions of interactions in powers of low-energy scales. In addition, they form a field- 
theoretic framework for extending the older effective-range expansions developed by Bethe 
and others [4| . In the context of charged particle scattering, Kong and Ravndal showed 
how Bethe's Coulomb-modified effective-range expansion could be represented as a pionless 
EFT. 

The results of Kong and Ravndal at NLO showed unexpectedly strong dependence of 
the scattering length on the subtraction scale fi. More recently, this approach was extended 
to NNLO and similar dependence of other coefficients on /i were observed, along with 
dependencies on the renormalization scheme and the choice of momentum structures in the 
NNLO terms. These results contrast with ones obtained from an analysis using a Wilsonian 

n 

renormalization group [8| , which showed a one-to-one correspondence between terms in the 
effective potential and the effective-range expansion. 

The main difference between the effective short-range potentials used in those two ap- 

nriQ 

proaches is that those in Refs. [5U6|, |7|] depend on the off-shell relative momenta of the two 
particles whereas those in Ref. |8| depend on energy. Applications of the renormalization 
group (RG) [9] to pure short-range interactions [10|] have shown that, in systems with bound 
or virtual states close to threshold, momentum- and energy-dependent terms have different 
scale dependence and hence appear at different orders in the power counting. 

To examine whether this differences between these short-range potentials is responsible, 
we have made an RG analysis of their scaling properties. To facilitate the comparison with 
Refs. we have used dimensional regularization (DR) with the power divergence 

subtraction (PDS) UJ, and modified minimal subtraction (MS) schemes. We consider 
potentials that can depend on energy to all orders but, for simplicity, we restrict our anal- 
ysis to the simplest momentum-dependent term, which is of second order in the off-shell 
momenta. 

The paper is organized as follows. In Sec. II, we apply the RG to an energy-dependent 
short-range interaction in the presence of the Coulomb potential. We concentrate on the 
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expansion around a marginally unstable nontrivial solution that is the analogue of the fixed 
point discussed in Ref. [10|. This allows us to confirm that DR leads to the same RG fiow 
as found in Ref. jsj using a sharp cut-off on the distorted waves. Then, in Sec. Ill, we find 
the leading momentum-dependent perturbation around our nontrivial solution. This has a 
different RG eigenvalue from the corresponding energy- dependent term. Its form suggests 
that it contributes only to the off-shell dependence of the scattering matrix, and we conform 
this by calculating the on-shell scattering amplitude in Sec. IV. Finally, we discuss the 
implications of this and other features of our results in Sec. V. 



II. SCATTERING BY LONG- AND SHORT-RANGE FORCES 

The two-body potential V between two charged particles can be separated into two pieces: 

V = Vc + Vs, (1) 

where Vc is the long-range Coulomb potential and Vs is the short-range strong interaction. 
For momenta that are too low to resolve the range of the short-range interaction, we can 
represent this part of the potential by delta functions and their derivatives. The RG will 
then lead us to a systematic power counting for these terms. 

For the sum of the two potentials above, the scattering matrix T can be split into Tc, 
the pure Coulomb T matrix, and Tsc which describes scattering between distorted waves of 
the Coulomb potential. If we restrict our attention to S-waves, the on-shell matrix element 
of Tsc can be related to the strong- interaction phase shift, Ssc{p), by 

Att p2i(To{p) 

where (Tq is the Coulomb phase shift, \ipcip)) i^" out-going Coulomb wave 

functions, and p = \/ME is the on-shell momentum.^ 

At low energies, the 5- wave phase shift 5sc(p) can be written in the form of an effective- 
range expansion as 

C{k/p)p cot 6sc{p) + 2nRe{H{K/p)} = -— + Kop^ + ■■■ , (3) 

ac / 



^ Here we are considering scattering of two identical particles of mass M. More generally M should be 
replaced by twice the reduced mass. 
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where 

27in 1 
= ^5^^ ' ^(^) = ^(^^) + ^ - ' 

K is the inverse of the Bohr radius (we will give its definition below), and is the logarithmic 

derivative of the F function. In this expansion, ac is the scattering length and ro is the 

effective range. 

This T matrix satisfies the Lippmann-Schwinger (LS) equation 

Tsc{E) = Vs + Vs GciE) Tsc{E) . (5) 

Here Gc = 1 / {E — Hq — Vc + i^) is the Coulomb Green's function. In terms of the Coulomb 
wave functions, this has the form 

J [2txY — + le 
III. RG ANALYSIS FOR ENERGY-DEPENDENT POTENTIALS 

Before examining momentum-dependent short-range potentials, we first consider purely 
energy-dependent ones. ™= ana,... .produce, t.e of Ref. g B.. w.M„ 

framework of DR and the PDS scheme. This establishes the notation we shall use for more 
general potentials and it allows us to make direct contact with the approach of Kong and 
Ravndal [s,^], which was also used in Ref. [3]. 

In DR, we replace the 3 + 1-dimensional Coulomb Green's function by 

Gc{E, /i) - M J p2_q2^,, ■ (7) 

The resulting integrals have poles at = 3 and 4, corresponding to linear and logarithmic 
divergences respectively. We subtract all of these, using the MS scheme for the poles at 
d = A. 

For a simple 5-function interaction, its strength Vs{p, fi) can depend on three scales: 
the on-shell momentum p, the inverse of the Bohr radius, 

aM 

«=^, (8) 
and the subtraction scale fi. The resulting regularized LS equation for Tsc takes the form 

Tsc{p, n) = Vs{p, n) + Vs{p, K, n)J^^{p, n)Tsc{p, k) , (9) 



where the bubble integral is defined by 



■^0= / 7^7T^(k|Gc;|k'). (10) 



(27r)3 (27r) 

In the dimensional regularization scheme used here, this becomes 



M kM 



where 7 is the Euler-Mascheroni constant. Note that the Coulomb wave functions at the 
origin, J (i^k/(27r)^(k|'?/'c'(g)), are finite and so the only divergences here arise from the 
integral over q. With momentum-dependent interactions, we encounter the derivatives of 
these functions. These diverge at the origin and so require additional regularization, as 
discussed in the next section. 

We now demand that the T matrix be independent of the arbitrary scale fi. This condition 
ensures that physical scattering observables will not depend on yU. It leads to a differential 
equation for the effective potential, which is similar to that controlling the cut-off dependence 



of the potential in Wilsonian approaches [8|, Il0|. This equation has the form 

^ Vsip, K, ^) = -Vs{p, K, ^ J^ip, ^^) = -Vsip, ^^Y (- ^ + • (12) 

This equation can be converted into an RG equation by expressing all dimensioned quan- 
tities in units of /i. We define p = p/yU and k = nj \i and introduce the rescaled potential 

ys{P-, /i) = ^ Vsifxp, fxk, /i) . (13) 

In terms of these, we can rewrite Eq. ( IT2l) as 

dfx \VsJ dp \VsJ dn \Vs J Vs 
Below all production thresholds, the effective potential should be an analytic function of 
the energy and all other low-energy scales. We thus require solutions to this equation to be 
analytic in p"^ and k. 

Fixed points of the RG equation play a crucial role, as discussed in detail in Ref. 10|. 
In the limit /x — 0, the rescaled potential, being dimensionless, is expected to become 
independent of fi. In other words, Vs should flow towards an infrared flxed point as yU — > 0. 

The only true flxed-point solution to Eq. (fT^ is the trivial one, Vq = 0. Although there is 
no other solution that is independent of n and analytic in k, there is a marginally unstable 



solution that flows logarithmically towards the trivial point. This behavior is the same as 
was found in [8] using a sharp cut-off on the basis of distorted waves. In the current scheme, 
it has the form 

-r^ = -l+2k\n(^) , (15) 

where A is another arbitrary scale. 

We now consider perturbations around this potential that scale with definite powers of 
H. Taking a single term of this form, 

=^-^-C/iV(p,^), (16) 

and inserting it into the RG equation (|T^ . we find that f{p,k) satisfies the eigenvalue 
equation 

Solutions to this are just products of powers of the low-energy scales, 

fip,k)=p'"^k\ (18) 

and their eigenvalues are 

u = 2m + n- l. (19) 

These include one unstable perturbation with a negative eigenvalue (m = n = 0), just as in 
the expansion around the nontrivial fixed point of the pure short-range potential [loj . There 
is also a marginal perturbation with z/ = (m = 0, n = 1) as expected from the logarithmic 
dependence on fi noted above. In the corresponding power counting, these terms can be 
assigned orders d = u — 1 = 2m + n — 2. 

The full short-range potential, expanded around the nontrivial solution Vq, is given by 

1 _ 1 



J2 C^nO/^'"^+""y"^/i" . (20) 



The coefficient of the marginal perturbation, Coio(A) depends logarithmically on A so that 
the full potential is independent of that arbitrary scale. 

To interpret the terms in this potential, we solve the LS equation with it and calculate 
the on-shell T matrix. Returning to physical units, the potential becomes 

^ ^ E ^-op'"^«" , (21) 



where 



1 



M 
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The corresponding scattering amphtude is 

{^ciP)\Tsc\i^ciP)) = 



— /i + 2k In 



C{k/p) 



(22) 



(23) 



Vs{p, K,/i) 



This leads to 

C{k/p) p cot 6sc + 2KRe{H{K/p)} = ^ C^^noP^™/?" + 2k 

m,n>0 

since the imaginary part of H is 

lm{H{r])} 



(24) 



C(r^) 



2r] 



(25) 



Comparing this with the Coulomb effective-range expansion, Eq. ([3]), we see that the 
coefficients of the energy-dependent perturbations are directly related the terms in that 
expansion by 



1 

ac 

1 



Cooo + Coio{A)k + 2k 
Cioo + Cno^ + 0{k^) 



f - 7 + In 



A 
4^ 



0{k' 



(26) 



If the coefficients in the potential are fine-tuned such that I/Vq exactly cancels the real part 
of then the T matrix develops a pole aX p = 0, corresponding to a zero-energy bound 
state. This is the scale-free system, analogous to the one corresponding to the fixed point 
in Ref. lOj. In the present case it is unstable against introducing a nonzero value for Cqoo 
or changing Cqio from the value that gives a vanishing coefficient for k. 

The full expansion around this point is a double one in powers of the energy (or p"^) and 
a (or k). For two protons k is approximately 3 MeV. This is small compared to and so 
this power counting can be used in the context of a pionless effective theory. More generally 
this scale is given by k = ZiZ2aMr, where is the reduced mass and the expansion in 
powers of k can break down if it is comparable to the scales of the underlying physics. For 
example, k for the system of two a particles is about 60 MeV, and so is of the order of m^r. 
In such be treated as a high-energy scale and the resulting pionless effective 



theory can be expanded in powers of energy only, as pointed out by Higa et al. 
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IV. RG ANALYSIS FOR MOMENTUM-DEPENDENT POTENTIAL 



We now turn to more general short-distance potentials that depend on momenta as well 
as energy. For simplicity, we consider potentials of the form 

Vs{p, K, k', k, /i) = Vi{p, K, n) + V2{p, K, n)k'^ , (27) 

where Vi and V2 are energy-dependent functions and k is the initial relative momentum. 
Although this expression is not Hermitian, it can easily be made so by adding a matching 
term with k — >■ k' . 

If we write the corresponding (off-shell) T matrix as 

Tsc = Ti + , (28) 

the (regularized) LS equation becomes two coupled linear equations 

= Vi + Vijf^T^ + V2J?^T^ , 
T2 = V2 + Vijf^T2 + V2jf^T2 . (29) 

Here Jq is the bubble integral defined above and 

J2= ! 7^7^A;2(k|Gc|k'). (30) 



(27r)3 (27r)3 

We apply DR and PDS to the integrals over q in the same way for both Jq and J2. 
However J2 contains an additional divergence, in the integral over k. This is because the 
Coulomb potential leads to wave functions whose derivatives are singular at the origin. 
We use DR and PDS for this too but with a different scale, A. Using the same scale to 
regulate all the integrals leads to an RG equation that contains nonanalytic dependence 
on the low-energy scales and is not suitable for a scaling analysis. The second scale A can 
be regarded as analogous to the factorization scale used in parton distributions to separate 
off the nonperturbative regime of QCD. Although the factorization scale can, and often is, 
taken to be equal to the renormalization scale, this is not a requirement. The Wilsonian 
approach to scattering in Ref. [sl used a 5-shell form for the short-range potential to avoid 
the nonperturbative singularities of several long-range potentials. The resulting physical 
amplitudes did not depend on the radius chosen for this. In the present case, results should 
be independent of the new scale A, and this will provide an important consistency check on 
our treatment. 
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With this regularization, we have 



Jfs ^ _ _ 2^A) Jo""' - AJ2 , (31) 
where is given by Eq. (fTTj) and 

AJ2 = ^/^V + 27rMft:3C'(-2) , (32) 

and the derivative of the Riemann zeta function is ('{—2) = —0.0304 • • ■ . 

By taking appropriate hnear combinations of the coupled equations for T12, we can 
rewrite them as equations that each involve only one piece of the potential,^ 



T2 = ¥2(1 + J^Ti + JpTa) . (33) 



Demanding that the T matrix be independent of fi, dTi 2/dfi=0, now leads to 



= ^{^ + J^'T, + J^'T2)+V,^T, + V,^T2, 

= ^{^ + j¥^T, + jf^T2)+V2^T, + V2^T2. (34) 

Multiplying these from the right by the factor (1 + J^^Ti + J^^T2)^^ and using Eq. ( 133|1 
for T1.2, we arrive at two coupled differential equations for Vi^2'- 

dV^ = _V^^_V,V2^, 

We rescale the potential as before, defining 

Viip, k, /i) = ^ Viiup, fik, iJ,) , V2ip, k,iJ,) = V2{fip, /i/l, fi) . (36) 
The rescaled functions Vi^2 then satisfy the coupled RG equations 

/x^vi = p—v, + k—v, + v, + {i-2k)v,' 



^ These are also easily obtained from the LS equation in the form Tsc — Vs + Tsc Gc Vs , instead of Eq. (O . 
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+ 



(1 - 2k) - 2k'^ - 2k 



K 



V1V2 



(37) 



r\ f) 

f^TrV2 = PTr^V2 + k—V2 + 3V2 + {i-2k)v,V2 

OjJ, op OK 



+ 



^2 ■ 



(38) 



We are interested in perturbations around the nontrivial solution described in the previous 
section, so we expand the potential around Vq as 



Vi = Vo + 5V, 



(39) 



and keep terms to first order in SVi and V2- The resulting linearized RG equations are 



OjJ, 



Op OK 



+ 



/i 



6Vi 

-2 



[l^2k){p^-2k^-2k-] +^k^ 



V0V2 



2 , 



Pgp+^Q^ + 3 + (1-2k)K) 



Vo 



(40) 
(41) 



The second of these is a homogeneous equation for V2 which does not involve 5Vi. By 
comparing it with the equation satisfied by Vq, we find that it has a solution 



V2 = CooiyuVo 



(42) 



with RG eigenvalue u = 2. Other solutions with larger eigenvalues can be obtained by 
multiplying this by powers of p"^ and k. 

Inserting this solution for V2 into Eq. (140!) . we get an inhomogeneous equation for 6Vi. 
This has the solution 



5Vi = -Cooi 



f -2k:^ -2k- ]Vo- — k^V^ 



IT 

~3 



+ fi^ {Ak^ + Bkp") . (43) 



The final term here is a solution to the homogeneous part of the equation, with the same 
RG eigenvalue as ¥2- It corresponds to two of the terms in Eq. (1201) when that is expanded 
to first in order in deviations from Vq. The full momentum-dependent perturbation has the 
form 



5Vs = Cooi 



P -p'^ + 2k^ + 2k-\Vo + — k^V^ 
Hj 3 



(44) 
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As already noted a Hermitian potential can be formed by adding the same structure with k 
replaced by k'. 

It is worth commenting on several aspects of this term in the potential. First, having 
solved the RG equation for fixed "factorization scale" A, we are now free to choose X = fi so 
that the whole expression is proportional to /i^. The RG eigenvalue is u = 2 and so the term 
is of higher order than the corresponding energy-dependent one, Cioofip"^ (which has u = 1). 



This is the same pattern as was observed in Ref. [10| for perturbations around the nontrivial 
fixed point of a pure short-range potential. As in that system, the momentum-dependent 
perturbation contains the structure k'^ — p'^ plus terms arising from the potential. This form 
suggests that it will vanish when acting on on-shell wave functions and hence it will alter 
only the off-shell behavior of the T matrix but not observables. We shall see below that this 
is indeed the case. 



V. SCATTERING BY MOMENTUM-DEPENDENT POTENTIAL 

We now take the momentum-dependent perturbation constructed in the previous section 
and combine it with the energy-dependent potential of Eq. fl2Tl) . which we now denote by 
VsE- Working in physical units, the potential has the form in Eq. ( |27|) . with 

Vi = VsE + SVu V2 = CooiVo, (45) 

where 

SVi = -Cool 



3 4:71 



+ {Ak^ + Bnp^) ^ V^f . (46) 



Using the LS equation in the form of Eq. (129|) . we can now calculate the T matrix. Expanded 
to first order in 6Vi and V2, this gives 

l/VsE-J^' [l/VsE-J^y \ViE VsE ' 
1 Vo 

T2 - (47) 

The on-shell T matrix can be written in terms of these as 

{i^c{p)\Tsc\i^t{p)) = ^0^1 + ^o^2T2 , (48) 
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where, following Refs. p, l6|, l7|, we have introduced 



(27r)3 
(27r)3 



A;2(k|^+(p)) = (/-2/s:2-2A/s:)^o. 



(49) 



Note that we have used DR to regularise the divergent derivative of the wave function in 
ip2i with the same scale A as used above for the similar divergence in J2. Inserting the 
expressions for Ti 2 from our potential gives 



C{kIv) 



{rc{v)\TscWc{v)) 



1 



^IVSE - Jf"" ■ V ^I^SE - } \ ^SE ySE WO 

where we have used the relation J2'^o ~ <^o^2 = —^J24'o- 

If we now substitute our explicit expressions for 6Vi and V2, we get 

„-2i(To 



5Vi V2 tp2 

VsE ^0 



1/2 
^SE 



+ 



v. 



AJ2 



SE 



(50) 



C{k/p) 



(V^^(p)|T5c|V^J(p)) 



SE 



+ 



SE 



SE 



12 " \VsE 
-2'kCooiMkW-2) 



(51) 



In solving the RG we kept only terms to first order in deviations from the fixed point. For 
consistency we should make the same approximation here, replacing Vss by Vq in the second 
term of this expression. The first term in the square bracket then vanishes. The remaining 
terms in cancel if we make the choices. 



A = 87r2CooiC'(-2) , B = 0, 



(52) 



for the constants of integration in our momentum- dependent perturbation. The on-shell 
scattering amplitude then reduces to the same form, Eq. fl2Tl) . as we obtained for the purely 
energy-dependent potential. This shows that the momentum-dependent eigenfunction of the 
RG equation is indeed an "equation of motion" term, affecting only the off-shell behaviour 
of the T matrix. 
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VI. DISCUSSION 



We have applied the RG to S'-wave scattering by a combination of Coulomb and short- 



range potentials, looking in particular at strong! 



close to threshold. Unlike previous work 



y interacting systems with bound states 
8j, our approach treats both energy- and 



momentum-dependence of the short-range potential separately. To handle the divergences 
generated by the contact interactions we use dimensional regularization with the PDS and 
MS schemes to subtract the linear and logarithmic divergences. The RG equation is obtained 
by requiring that the off-shell scattering matrix be independent of the subtraction scale. 

The resulting RG equation has a physically interesting nontrivial solution. This is not an 
exact fixed point because it possesses a marginally relevant perturbation and so it evolves 
logarithmically with the regulator scale. Perturbations around this can be expanded in pow- 
ers of energy, momenta and a (or the corresponding low-energy scale k = aM/2). The RG 
eigenvalues and hence the power counting for purely energy-dependent^erturbations are the 
same as those found for systems with only short-range interactions [10]. This reflects the 
fact that the Coulomb and free wave functions have the same power-law behavior near the 
origin. The fact that the Coulomb potential leads to a logarithmic discontinuity in deriva- 
tives of the wave functions is reflected by the presence of the marginal term proportional to 

K,. 

Like the fixed point of the pure short-range potential, the nontrivial solution describes a 
system with a bound state at zero energy. The double expansion around it in powers of the 
energy and k corresponds to the Coulomb effective-range expansion originally introduced by 
Bethe {4], with the coefficients expanded in powers of a. 

We have extended this approach to terms that depend on the off-shell momenta, present- 
ing in detail the analysis for the simplest of these, which is second-order in the momenta. 
The results for the power counting match with those found in the pure short-range case 
|lol | . (Again, this is not surprising given the behavior of the Coulomb wave functions at the 
origin.) In particular, the momentum-dependent terms appear at higher orders than the 
corresponding energy- dependent ones. The one studied here has an RG eigenvalue u = 2, 
compared with z/ = 1 for the term proportional to the on-shell energy. 

The structure of these momentum-dependent terms indicates that they arise from the 
equation of motion and hence they should affect only the off-shell behavior of the scattering 



13 



matrix. We have verified that our second-order term does not contribute to scattering 
observables. 

This RG analysis can explain some of the puzzling features seen in other applications 
of EFTs to charged-particle scattering. The different power counting for energy- and 
momentum-dependent perturbations around the nontrivial fixed point means that one can- 
not "use the equations of motion" to remove energy- dependence from the effective potential 
without mixing terms of different orders. Requiring the potential to be energy-independent 
for all values of the regulator scale is only possible if the higher-order momentum dependent 
perturbations are given unnaturally large coefficients. This is quite different from expansions 
around a trivial fixed point where the corresponding terms appear at the same order. 

Most approaches to the EFT for charged-particle scattering are based on purely 
momentum-dependent potentials (derivatives of 5-functions) . Each term in such a 

potential can be built out of a combination of energy-dependent perturbations and off-shell 
ones. The leading-order piece of each is the energy- dependent one and so these potentials 
can reproduce the Coulomb effective-range expansion. However, as mixtures of different RG 
eigenfunctions, these terms contain pieces that violate the power counting. As a result, the 
inclusion of higher-order terms can significantly change the coefficients of lower-order ones. 

Moreover, the same on-shell perturbation can form part of more than one distinct 
momentum-dependent term. For example, Ref. shows that both the structures k'^ + k"^ 
and fc^fc'^ can generate a shape parameter in the effective-range expansion. However they 
contain different admixtures of off-shell perturbations and hence they affect the renormal- 
ization of the lower-order terms differently, as demonstrated by the numerical results in that 
paper. 

Finally, the combination of momentum-dependent perturbations and the Coulomb po- 
tential leads to new divergences. These arise from derivatives of the Coulomb wave functions 
at the origin. We have used DR and PDS to regularize these but we found it necessary to 
keep their subtraction scale distinct from the one used for the integral over intermediate 
states. We regard this second scale as analogous to a factorization scale, separating off the 
nonperturbative regime caused by the singularity in the long-range interaction at the ori- 
gin. Although we have been unable to obtain a useful RG equation by demanding that the 
off-shell scattering matrix be independent of this scale, we are able to show that physical 
observables do not depend on it. 
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